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Vortex-induced disturbance field 
in a compressible shear layer 

By D. Papamoschou 1 AND S. K. Lele 2 


The disturbance field induced by a small isolated vortex in a compressible shear 
layer is studied using direct simulation in a converted frame. The convective Mach 
number, M c , is varied from 0.1 to 1.25. The vorticity perturbation is rapidly sheared 
by the mean velocity gradient. The resulting disturbance pressure field is observed 
to decrease both in magnitude and extent with increasing M c , becoming a narrow 
transverse zone for M c > 0.8. A similar trend is seen for the perturbation velocity 
magnitude and for the Reynolds shear stress. By varying the vortex size, we verified 
that the decrease in perturbation levels is due to the mean-flow Mach number an 
not the Mach number across the vortex. At high M c , the vortex still communicates 
with the edges of the shear layer, although communication in the mean-flow direc- 
tion is strongly inhibited. The growth rate of perturbation kinetic energy declines 
with M c primarily due to the reduction in shear stress. For M c > 0.6, the pres- 
sure dilatation also contributes to the decrease of growth rates. Calculation of the 
perturbation field induced by a vortex doublet revealed the same trends as in the 
single-vortex case, illustrating the insensitivity of the Mach-number effect to the 
specific form of initial conditions. 


1. Introduction 

It is well known that the Mach number has a powerful effect to suppress the in- 
stability and growth of free shear flows. Landau 1 first showed that the vortex sheet 
becomes stable when the relative Mach number exceeds a critical value (V2 for the 
equal-density case). Early linear analyses (Lin 1953, Gropengiesser 1970, Blumen e 
al 1975) and single-stream experiments (Sirieix k Solignac 1966) discovered tha 
the growth rates of compressible, finite-thickness shear layers decrease sharply with 
increasing Mach number. Recent linear analyses (Ragab k Wu 1988, Zhuang et 
al 1990), computations (Sandham k Reynolds 1989, Lele 1989), and two-stream 
experiments (Chinzei et al. 1986, Papamoschou k Roshko 1988, Samimy k Elliott 
1990 Goebel k Dutton 1991) covering a larger range of conditions confirmed the 
above trends. The experiments, in particular, showed that at high Mach numbers 
the turbulent shear-layer growth rate decreases to as little as one fifth of the in- 
compressible value. The perturbation levels of velocity decline in a similar fashion 
(Samimy k Elliott, Goebel k Dutton). Despite these recent gains m the held ot 
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compressible turbulence, the fundamental physical reason for the stabilizing effect 
of Mach number remains elusive. 

While compressible shear flows typically contain density gradients, density effects 
alone are not responsible for the large reduction in growth rates. Brown & Roshko 
(1974) showed that the growth rate of the subsonic, variable-density shear layer 
changed only by about 50% when the density ratio was varied by a factor of 50. 
There is, therefore, a large effect associated with the Mach number itself. An 
exclusive property of Mach number is the ability to cut off communication between 
parts of the flow, a well-known phenomenon in supersonic flow. 

Morkovin (1987) stresses that upstream and cross-flow communication is essential 
for instabilities at supersonic and hypersonic speeds. His point is primarily based 
on Mack’s (1984) linear stability analysis, where it is shown that the most unstable 
waves are those whose phase speed is subsonic relative to the free-stream velocity. 
Morkovin suggests the existence of zones of influence, defined by Mach cones, out- 
side which a disturbance is not felt. The concept of reduced communication at high 
Mach number has been incorporated into recent turbulence models. BreidenthaPs 
(1990) sonic-eddy model is based on the assumption that turbulent eddies whose 
rotational Mach number is greater than unity do not participate in fluid entrain- 
ment, while those with rotational Mach number of unity or less engulf fluid like 
incompressible eddies. The mixing-length model of Kim (1991) for a supersonic 
shear layer assumes that disturbances do not penetrate outside a region bounded 
by relative sonic velocities. 

Morkovin s concept and the models by Breidenthal and by Kim, while conceptu- 
ally useful, do not explain the fact that stabilization starts at low subsonic values of 
the Mach number and is practically complete when the velocity difference between 
the center and the edge of the shear layer becomes sonic (see references above). 
Furthermore, the fact that the relative velocity is supersonic does not prove that 
information will not propagate outside the region bounded by the sonic velocity. 
To address the problem in a more quantitative, although still idealized fashion, Pa- 
pamoschou (1991) used ray theory to study the acoustic field of a monopole placed 
inside a shear layer. It was shown that the influence of the monopole diminishes 
with increasing Mach number but that a low level of communication still exists 
at supersonic relative speeds. Still, the connection between communication and 
instability remained a speculative one. 

In this study, we explore the link between acoustic field and disturbance mo- 
tion by direct numerical simulation of a simple instability problem: the impulse 
response of a shear layer to a localized vorticity perturbation and the change of 
that response with increasing Mach number. Specifically, we examine the unsteady 
process by which the shearing of a single vortex, placed in the center of the shear 
layer, produces velocity, pressure and density fluctuations in the surrounding flow 
field. Since vortical interactions among eddies are an essential ingredient of turbu- 
lence, our simplified problem may serve as a building block for understanding the 
effect of Mach number on the more complex vortical interactions in realistic shear 
flows. 
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2. Initial conditions , 

The computation is two dimensional, inviscid and temporal, with periodic ^ bound- 
ary conditions in the mean-flow direction and non-reflecting boundary conditions in 
the y-direction. At time t = 0, the flow field consists of the shear layer described by 
a hyperbolic-tangent velocity profile U(y), with velocity difference A U and vorticity 
thickness 6 „ = 0.1, and an Oseen vortex (Fig.l). 


A U/2 




L = 0.1 5* 
r=- 0.01 A UL 


Shear layer 


+ Oseen Vortex 


papamos/figl-ps [1] Figure 1. Initial conditions. 

The Oseen vortex is a solution to the flow with zero shear and has a tangential 

velocity _ 
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where T is the circulation, L is the vortex size, and a = 1.256431 is chosen such 
that v« = vo, mai at r = L. The corresponding pressure distribution, derived by the 
radial momentum equation under the assumption of homentropic flow, is 

( 2 ) 
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with Ei(x) the exponential-integral function and oo denoting the unperturbed con- 

dl The circulation of the Oseen vortex is fixed here at T = 0.01 AUL. Unless other- 
wise stated, the size of the vortex is L = 0.U u . The speed of sound a is constant 
throughout the flow field. The convective Mach number is 

AU 
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with AU the velocity difference across the layer. To increase Af c , the velocity 
difference remains constant and the speed of sound decreases. The M c range covered 
here is from 0.1 to 1.25. 

For finite shear, the Oseen vortex is no longer a solution to the flow field. As 
a result, the flow tries to adjust to the new condition by generating an acoustic 
wave which propagates in all directions and becomes distorted by the mean Mach 
number gradient. Behind the acoustic wave, a disturbance field is established whose 
pattern evolves slowly with time. The disturbance velocity field is u = (u, v) and 
the pressure perturbation is p. 

The computation is advanced in increments of acoustic time t whose non dimen- 
sional version is t a = ta/6 w . The corresponding shearing time is 

t, = tS = t a — = 2 t a M c 
a 

where 5 = AU/S w is the maximum shear. At fixed t 9 , a particle at the edge of the 
shear layer has traveled the same distance, hence the flow has been sheared by the 
same amount, regardless of the value of M c . Comparisons among the flow fields at 
various A/ C ’s will be made at fixed shearing time. 

3. Computational Details 

The computations were carried out in a domain of size 2 in the mean-flow direction 
(x) and of size 0.5 in the transverse direction (y). The disturbance vortex was placed 
at the center of the domain which was discretized by a non-uniform mesh in both 
directions. In the computations reported in this paper, the mesh contained 200 
points in x (with a maximum mesh stretching of 3.2) and 50 points in y (with 
a maximum mesh stretching of 3.2). Near the vortex the mesh size was 0.0026 in 
both x and y . Spatial derivatives were evaluated using the sixth-order compact finite 
differences (Lele 1989, 1992), and third-order compact storage Runge-Kutta scheme 
was used for time advancement. Periodic boundary conditions in x and “non- 
reflecting” boundary conditions in y were employed. To make the initial conditions 
compatible with the periodicity in x, the method of images was used (with two 
images taken outside each x- boundary). To suppress the spurious generation of 
26 error waves in regions of mesh stretching, compact filtering scheme designed to 
remove only the near 2 6 waves (Lele 1992) was applied to the fields being integrated 
every 20 time steps. It was verified that the applied filtering produced no significant 
change in the well-resolved physical disturbances. 

4* Results and Discussion 

We first present the time evolution of several quantities at M c = 0.4. Fig. 2 
shows the vorticity deformation due to the mean velocity gradient. At late times, 
vorticity is concentrated in a narrow, almost horizontal layer. Even though only the 
evolution at M c = 0.4 is presented, the vortex shearing seen in Fig. 2 is practically 
the same for the full range of Me s covered here. The corresponding evolution in 
kinematic Reynolds shear stress -uv is depicted in Fig. 3. At t s = 0, the round 
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Figure 2. Iso-contours of disturbance vorticity uL/AU for M c = 0.4. Contour 
levels: Minimum=-0.075; Maximum=-0.005; Increment=0.005. x-increment=0.02; 
y-increment=0.01. 

vortex has equal amounts of positive and negative shear stress. As time progresses, 
the negative part vanishes while the positive part occupies a larger part of the flow 
field. The time development of the divergence is seen in Fig. 4, where the expanding 
wavefront is evident. As mentioned earlier, the wavefront arises from the reaction 
of the flow to the initial conditions. Partial reflection of the wave front from the 
upper and lower boundaries is due to numerical error associated with the boundary 
conditions. 

We now present the Mach number effect on the disturbance field at fixed t 3 . Fig. 
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i.= 0 



t 1 r 





FIGURE 3. Iso-contours of Reynolds shear stress — uv/(AU) 2 for M c = 0.4. Con- 
tour levels :Maximum = 26E- 5 ; Mini mum=-26E- 5; Increment=3E-5. ar-increment=0.1; 
y-increment = 0 . 05 . 


5 depicts the magnitude of the pressure perturbation at four different Af c ’s. The 
extent and magnitude of the field decrease rapidly with increasing Mach number. 
At M c = 1.25, the pressure field is reduced to a narrow transverse zone. The same 
trend is seen for the disturbance-velocity magnitude, shown in Fig. 6, and for the 
shear stress, shown in Fig. 7. Comparison of Figs. 5 and 7 shows that the decline 
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Figure 4. Iso-contours of divergence V • u L/AU for M c = 0.4. Contour 
levels: Maximum=0.025; Minimum=-0.017; Increment=0.004. x-increment=0.2; 
y * increment =0 . 05 . 

in magnitude and extent of the p and -uv fields is very similar. 

The pressure pattern seen in Fig. 5 reveals some very noteworthy effects of 
compressibility on the vortex influence. Evidently, at high M c , the influence of 
the vortex does not propagate in x but stays confined within a narrow transverse 
region. Interestingly, the vortex still communicates to edges of the shear layer, even 
when the relative velocity is supersonic. This implies that the turbulence models 
mentioned in the Introduction, which assume that all interactions occur within the a 
layer bounded by sonic velocities (Breidenthal 1990, Kim 1991), are overly idealized. 
The lack of communication in the x-direction may explain the lack of vortex pairing 
and shear-layer roll up at high M c ■ 

We now examine the source term for the transverse component of the kinetic 
energy ±t> 2 , namely -vdp/dy. Note that generation of \v 2 does not involve the 
mean shear, thus is directly related to the pressure field. Fig 8 shows the source 
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1 A/ c = 0.80 




Figure 5. Iso-contours of pressure disturbance magnitude |p|/(/»ooA^) at t, = 
4.2. Contour levels: Maximum=0.032; Minimum=0.002; Increment— 0 .002 x- 
increment=0.2; y-increment=0.05. 


term dramatically decreasing with increasing M c , practically vanishing at M c = 0.8. 
At M e = 1.25, the radiative nature of the flow, coupled with boundary-condition 
errors, make this term reappear in roughly-equal negative and positive amounts, so 
its integrated contribution is near-zero. The decline of —vdp/dy, and the resulting 
reduction in |w|, are the direct result of the reorientation of the pressure field due 
to the Mach number. Since |u| and | v | are of the same order, the shear stress —uv 
is also reduced, as seen in Fig. 7. The lower —uv, in turn, causes a decrease in 
the growth rate of the overall fluctuation kinetic energy, which will be shown later. 
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FIGURE 6. Iso-contours of disturbance velocity magnitude |u|/AC/ at t, - 
4.2. Contour levels: Maximum=0.028; Minimum=0.002; Increment— 0.00-. x 
increment =0.2; y-increment=0.05. 


The same mechanism has been observed to reduce the kinetic-energy growth rate 
in compressible homogeneous sheared turbulence (Blaisdell et al. 1991). 

To verify that the observed trends are due to the mean-flow Mach number and 
not due to the Mach number M L = M C L/6 U across the vortex, we varied Ml by 
changing the vortex size, keeping the circulation fixed, and compared the resulting 
disturbance fields. Fig. 9 offers two such comparisons: at M c - 0.4, Ml changed 
from 0.04 to 0.08, while at M c = 0.8, Ml changed from 0.08 to 0.04. If Ml were 
the critical parameter, one would expect similarity between the fields at Ml - 0.04 
and between the fields at Ml = 0.08. This is clearly not the case, and, even though 
details near the vortex core change with size, the overall extent and magnitude o 


the fields depend only on M c . .... 

A quantity of great interest in any instability problem is the fluctuation kinetic 
energy fc = ±u • u and its generation terms. For our inviscid, two-dimensional flow, 
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Figure 7. Iso-contours of Reynolds shear stress - uv/(AU ) 2 at t a = 4.2. Contour 
levels: Maximum=33E-5; Minimum=-5E-5; Increment=2E-5. Negative values are 
denoted by dotted lines. x-increment=0.2; y-increment=0.05. 

the instantaneous k is described by 
dk dk dU 

dt +U di = ~ uv 'd^ +pVu - v - M - v • ( ku ) (3) 

The terms on the right-hand side are called production due to shear stress, pressure 
dilatation, pressure transport, and kinetic-energy transport. We examine the time 

evolution of the quantities in Eq. 3, integrated in space over our computational 
domain. 
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FIGURE 8. Iso-contours of (vdp/dy)/ Poo AU 3 at t, = 4.2. Contour levels: Maxi- 
mum=0.0026; Minimum=-0.0022; Increment=0.0002. Negative values are denoted 
by dotted lines. i-increment=0.2; y-increment=0.05. 

The growth of the disturbance kinetic energy versus time for different M c ' s is 
seen in Fig. 10(a). As expected, the higher M c ’s produce significantly lower growth 
rates. At early times, the growth rate of the M c = 1.25 case slightly exceeds that of 
the M c = 0.8 case. This is probably due to the longer presence of the wavefront in 
the computational domain at M c = 1.25. Recall that M c is increased by lowering 
the speed of sound, thus the wavefront and its velocity perturbation take longer time 
to propagate outwards at high M c . By late times, however, the M c = 1.25 growth 
rate is seen to saturate and decline below the M c = 0.8 curve. Figure 10(b) shows 
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M c = 0.80 
!/«„ = 0.1 




Figure 9. Iso-contours of Reynolds shear stress - uv/(AU ) 2 at t, = 4.2 for 
variable initial vortex size L. Contour levels: Maximum=26E-5; Minimum=-26E-5; 
Increment=3E-5. Negative values are denoted by dotted lines. x-increment=0.1; 
y-increment=0.05. 


the production due to shear stress, which also declines rapidly with increasing AI C . 
The pressure dilatation term, seen in Fig 10(c), is much lower than the production 
term for M c = 0.1 and 0.4, but becomes comparable at the higher M c 's. Thus, the 
decline in growth rates, most of which occurs at subsonic M c , is primarily due to 
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FIGURE 10. Evolution of spatially-integrated quantities versus t, for variable M c : 
(a) perturbation kinetic energy; (b) production due to shear stress; (c) pressure 
dilatation, o M c = 0.1; A M c = 0.4; + M c = 0.8; x M c = 1.25. 


the reduction in shear stress. The pressure dilatation term becomes significant at 
high-subsonic and supersonic M c and makes the overall generation of kinetic energy 
even smaller. The pressure transport and kinetic-energy transport terms were found 
to be negligible compared to the previous two terms. 

To test the sensitivity of the Mach-number effect on initial conditions, we replaced 
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Figure 11 . Iso-contours of dipole disturbance vorticity ljL/AU for M c = 0.4. 
Contour levels: Minimum=-0.07; Maximum=0.07; Increment=0.01. Negative val- 
ues are denoted by dotted lines. x-increment=0.02; t/-increment=0.01. 


the single vortex by a vortex dipole. The dipole vortices have the same circulation 
and size as the previous single vortex, and are separated by a distance 4 L. The 
shearing of the dipole versus 1 3 is shown in Fig. 11. At late times, the positive 
vorticity is stretched a little more than the negative one, probably because of the 
interaction between the two vorticity regions. The Reynolds shear stress versus M c 
is depicted in Fig. 12, where its magnitude and extent are seen to shrink rapidly with 
increasing M c , as in the single- vortex case. Additionally, we reversed the circulation 
of the single vortex, producing a counter-clockwise motion. The magnitude and 
extent of perturbation fields remained identical to those of the clockwise vortex. 
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FIGURE 12. Iso-contours of dipole Reynolds shear stress -uv/(AU) at t, - 4.2. 
Contour levels: Maximum=21E-5; Minimum=-5E-5; Increment -2E- -5. Nega ive 
values are denoted by dotted lines. x-increment=0.1;y-mcrement-0.05. 


Evidently, the precise form of the initial condition does not affect the stabilizing 
influence of the mean-flow Mach number. 

To illustrate the robustness of the high-M, features observed here, we make qual- 
itative comparisons among pressure fields at M, = 1.25 generated by three Afferent 
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Figure 13. Qualitative comparison of M c = 1.25 pressure disturbance fields Di- 
rect simulation:(a) vortex monopole; (b) vortex dipole; (c) acoustic monopole. Geo- 
metric acoustics theory:(d) acoustic monopole. ar-increment=0.2; j/-increment=0.05. 

initial conditions. In Fig. 13, the pressure-disturbance fields for (a) a single vor- 
tex, (b) a vortex dipole, and (c) an acoustic monopole are displayed. The acoustic 
monopole has an initial-pressure field same as in the single Oseen vortex, but the 
vorticity is now zero. All three features look very similar, which once again shows 
heir dominance by the mean -flow Mach number and their independence from de- 
tails of the initial condition. Fig. 13(d) depicts the field of an acoustic monopole, 
computed by the simple energy- invariance method of geometric acoustics (for details 
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see Ref. 19). Even though geometric-acoustics theory is only applicable to high- 
frequency sources, the calculated pressure field resembles very much those computed 
by direct numerical simulation. 

5. Conclusions 

The disturbance field generated by the shearing of a vortex placed in the cen- 
ter of a compressible shear layer has been studied by means of direct numerical 
simulation. The magnitude and extent of the induced pressure and velocity fields 
diminish rapidly with increasing convective Mach number M c , with most of the de- 
crease occurring at subsonic values of M c . At M c > 0.8, the influence of the vortex 
is reduced to a zone of narrow streamwise extent. The vortex still communicates to 
the edges of the shear layer, although communication in the mean-flow direction is 
practically cut off. Perturbation kinetic-energy growth rates decrease with increas- 
ing M c , primarily due to reduction of the Reynolds shear stress. Pressure dilatation 
contributes to the decline in growth rates for M c > 0.6. The trends observed here 
are insensitive to the specific form of the initial disturbance and are a result of the 
mean-flow Mach number. 
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